Low-temperature transport in Heisenberg chains 



J.V. Alvarez and Claudius Gros 
Fakultat 7, Theoretische Physik, University of the Saarland, 66041 Saarbriicken, Germany. 



(N 
O 
O 
(N 

C 

a 

Hi 

OO 



> 
in 
00 
in 
in 
o 



c3 



o 
o 



A technique to determine accurately transport properties 
of integrable and non-integrable quantum-spin chains at finite 
temperatures by Quantum Monte-Carlo is presented. The re- 
duction of the Drude weight by interactions in the integrable 
gapless regime is evaluated. Evidence for the absence of a 
Drude weight in the gapless regime of a non-integrable sys- 
tem with longer-ranged interactions is presented. We estimate 
the effect of the non-integrability on the transport proper- 
ties and compare with recent experiments on one-dimensional 
quantum-spin chains. 

PACS numbers: 75.30.Gw, 75.10.Jm, 78.30.-j 



Introduction - During the last few years several fam- 
ilies of materials containing well characterized quasi one- 
dimensional spin-1/2 structures have been synthesized. 
The charge-transfer gap is in many cases large and the 
spin excitations contribute significantly to the thermal 
and magnetization transport at low temperatures. For 
example, 63 Cu NMR studiesEl in S^CuOs have mea- 
sured a spin diffusion coefficient (equivalent to diffu- 
sive magnetization transport) several orders of magni- 
tude larger than the value for conventional diffusive sys- 
tems, and thermal transport measurements in Sr2Cu03 
and SrCu02 indicates quasi-ballistic transport with a 
mean-free path of several thousands of A. 

These unusual results have been related to the peculiar 
physics of one-dimensional quantum chains, ft is known 
that the spin transport in the XXZ chain 
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A connection be- 
enlights the situa- 



is not diffusive, even for T 3> J x , 
tween integrability and transport! 

tion. In a generic integrable model, like H^ xxz \ the spin 
current is not conserved, but it has a non- vanishing 'com- 
ponent' with respect to the projection onto a conserved 
operator. As a consequence, the current-current correla- 
tion functions do not decay to zero for large times. 

The hamiltonian appropriate for real compounds like 
Sr2Cu03 and SrCu02 correspond to H^ xxz ' only in first 
approximation. A crucial question is therefore: How 
do small deviations from integrability change the pic- 
ture described-_above? Extending an earlier analysis 
by GiamarchiEHI, Rosch and Andreij-concluded recently 
within a memory-matrix approachtil, that deviations 
from integrability lead to an exponentially large conduc- 
tivity in Hubbard-like models away from commensurabil- 
ity 

Despite the on-going effort devoted to this problem, 
the fundamental difference in between integrable and 



non-integr able models have not yet shown up in QMC- 
simulationai3il3, presumable due to the demanding nu- 
merical requirements. In this Letter we develop the tech- 
niques and tools for the data-analysis necessary for ex- 
tracting the Drude weight for integrable and the life- 
time, mean-free path and the diffusion constant for non- 
integrable quantum-spin models from correlation func- 
tions obtained by QMC simulations using the Loop 
AlgorithmEa. We present extensive comparison with 
exact results obtained by Bethe-Ansatz and compare 
for a non-integrable system with recent experiments on 
Sr 2 Cu0 3 . 

Drude weight - QMC-simulations yield in general 
correlations functions on the imaginary-time axis. We 
therefore consider the Kubo formula for the dynamical 
conductivity a(u> n ) — lim g _>o c(<7, w„) 



- (K) - A(q,u> n ) _ D(q,w n ) 



(1) 



on the imaginary axis, where (K) , for the models we will 
consider, is the expectation value of the kinetic energy 
per site and where A is the current-current correlation as 
a function of the Matsubara frequency, 



A(q,u r , 
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(f(q,r)f(-q,0))dr. (2) 



Eq. (0) leads via a(ui) = t:D{T)6(lu) + a reg (cj) to the 
conventional relationEj between the current-current cor- 
relation function and the Drude weight (we are bor- 
rowing the terminology of electrical responses) D(T) = 



lim w „_ lim g _ D(q, lo Ti 



D(T) = -(K) - A(g^0,w n ^0) . 



(3) 



Note that A(q, 10 + id) is analytic in the upper half of the 
complex w-plane and the extrapolation along the imagi- 
nary axis can be reliably performed at low temperatures, 
when many Matsubara-frequencies u> n —Jh:Tn are avail- 
able close to ui = for the extrapolationEj. 



The continuity equation -^S[(t) 



to the expression 



&jf (t) = leads 



1+1 



(4) 



for the current operator ji(t) in the XXZ-model. In 
Fourier space, the continuity equation takes the form 



S*(t) = [H,S'] = i(l 



(5) 



A(g, u> n ) is a non-diagonal four-site operator. In principle 
non-diagonal operators can be computed using the loop 
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algorithnJl! but, as it was discussed recentlyQ, the algo- 
rithm to compute two-site correlation functions is more 
efficient than the algorithms that compute operators in- 
volving three or more sites. To obtain high quality data 
we measure directly the dynamical susceptibility in imag- 
inary time 



S(q,LO r , 







> t (S'Jt)S'(P)) , 



(6) 



a simple diagonal two-site correlation function. The 
Drude peak is computed by using Eq. ([^) and the re- 
lation 



■>nS{q,(- 



Asm 2 (q/2)D(q,u ri 



(7) 



This relation can be derived by a two-fold partial inte- 
gration of the right-hand side of Eq. (||) with respect to 
r, which leads to 

S(q,u n ) = ^([[H,S*],S*_ q ]) - 4sm ^ g/2) A(g,o; w ) , 



where we have used Eq. (||) for the imaginary-time deriva- 
tives of Sq and St q and definition (^) . Evaluation of the 
double commutator [[H, S q ],St q ], which is the boundary 
term from the partial integration, then leads to Eq. (0). 
Note, that this double commutator occurs here for the 
Matsubara correlation functions and does not occur for 
a related real-frequency correlation functional. 

Data- Analysis - At low temperatures and frequen- 
cies, the scaling of D(q,uj n ) can be obtained simply in- 
voking the conformal symmetry of the model emerging 
in the gapless regime J z < J xx . S(q,uj n ) at small q takes 
then form S(q } u! n ) = D x (T)q 2 /((cq) 2 +u%). This expres- 
sion and Eq. (Q) suggest the form 



D(q,uj r , 



A 2 (g)+u; 



(8) 



Alternatively, Eq. (||) can be viewed as the first term 
of the exact representation for D(q,w n ) containing an 
infinite-number of termscj. 

The XXZ-model maps to an interacting ID spinlcss 
fermionic system at half filling. For the noninteracting 
case (the XX chain) we can compute exactly D(q, oj n ) and 
we obtain D(0) = J xx Jir, and A(q) = J xx sxn(q) ~ J xx q. 
For J z < J xx the Umklapp term remains marginally ir- 
relevant ami one expects Luttinger Liquid like correlation 
functions^ like Eq. (§). 

In the ideal Luttinger Liquid the bosonic excitations 
are arbitrarily well defined at low temperatures and low 
frequencies. The conductivity is necessarily infinite in 
that case, and ansatz (||), being invariant under time re- 
versal reflects that property. Nevertheless, we are also in- 
terested in the study of more general situations in which 
the bosons can decay, become quasiparticles, and mem- 
ory effects can be taken into account. To this end we use 
the formalism of the memory matrices which has been 



successfully applied to study how jipne-dimensional elec- 
tron liquids can gain a resistivityEJO. 

This is achieved by taking the following fitting function 
for D(q,uj n > 0): 



D(q,v n ) = ^2 



f^A 2 (q) + 2 lj (q)u n + 



-^7 



(9) 



D(q,u n ) is analytic in the upper complex-plane for 
lj{q) > 0. For Aj(0) = lim g ^o Aj(g) we find either a 
(i) gapless regime: Ai(0) = and A 2 (0) > 0, or a 
(i) gap-full regime: Ai(0) > and A 2 (0) > A x (0). 
The first term in Eq. ([)]) dominates the low-frequency 
behavior in both cases and we have set generally 72 = 
in order to keep the number of parameters to a minimum. 

In the gapless regime the optical conductivity ([|) takes 
for small frequencies the Drude form 



Recr(w) 



2Di(0)7i(0) 



+ 4 7l 2 (0) " l + {uOTf 



(10) 



where we introduced the DC-conductivity <7q = 
Z?i(0)/ (271(0)) and the quasi-particle lifetime r = 
(2 7l (0))" 1 . For r -» 00 Eq. @ reduces to Recr(w) = 
7rZ?i(0) S(u>). At high frequencies 

lim D(0,w n ) = -(K) = D 1 {0) + D 2 (0) (11) 

uj n — >oo 

and a finite /^(O) results in a reduction of the Drude 
weight D(T) with respect to the kinetic energy, see Eq. 
(^). A finite Z?2(0) measures therefore the amount of 
decay experienced by the total current due to the inter- 
actions. We note that the Ansatz Eq. ([)]) for D(q,u>rX, 
together with Eq. (pd|), is consistent with the f-sum-ruleE-D 
J °° Kea(uj) = — f (K) for the optical conductivity. 

Ballistic transport - In Fig. [I] we show the values for 
Di(q) and Ai(g) for theJCX chain for different system 
sizes as obtained by QMQL3, D2 (q) and 71 (q) optimize to 
zero for this model. The different system sizes collapse 
at the smallest momenta and in this way the thermody- 
namic limit and the q — > limit are performed simultane- 
ously. Our prescription to extract -Di(O) as the limiting 
value D\ (0) = lim g _,o Efeiq) is more involved than pre- 
vious q = calculations^! but we believe it to be nec- 
essary in order to avoid substantial finite-size effects due 
to the fact that D(q — > 0, uj n ) becomes a step-function at 
u> n — in the thermodynamic limit, compare Eq. (||). 

For finite interaction J z we find in the gapless regime 
Jz < Jxx, that Ai (q) ~ c( J z ) q where c( J z ) is the velocity 



c(J z ) 



n yJxx Jz 



(12) 



2 arcos( J z /J xx ) 

of the des Cloiseaux-Pearson spectrumEi. In Fig. || we 
compare the measured value of c{J z ) with the Bethe 
Ansatz result @ for J z /J xy = 1/2,1. A x {q) is fitted 
well in the gapped phase by e(q) — yAo + ( c l) 2 - We 
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find Ao = 0.191 J xx , widen is close to twice the one- 
magnon gap of 0.091 JxJlj. 

The damping 71 (q) is vanishing small for J z < J xx and 
acquires a finite value in the gapped phase which can be 
fitted phenomenological by the relation 7i(q)Ai(g) « 
const., independent of q. 

In Fig. H (inset) we present values obtained by QMC 
for the Drude weight in the gapless regime, J z < J xx at 
T = 0.004 J xx . We find good agreement with the T = 
Bethe Ansatz resulted 



D(0) = 



v 2 smiixjv) 
47T v - 1 



Jz = Jxx cos(-k/v) , (13) 



specially at the smaller values of J z . In this region the 
correlation function is expected to satisfy more accu- 
rately the scaling law (||). The S{q,io n ) in the isotropic 
Heisenberg model does not follow strictly the Luttinger 
Liquid form. It presents multiplicative logarithmic cor- 
rections, that should be observed also for values of J z 
slightly smaller than J xx when finite temperatures and 
system sizes impose a cut-ofLfor the RG equation of the 
Umklapp coupling constant^]. 

We study now the behavior of the Drude weight at 
finite temperatures for models free from strong multi- 
plicative corrections. The main conclusion of Zotoscj is 
the fast decay of the Drude weight when the tempera- 
ture increases, D(T) — D(0) ~ —kT 17 ^ where n is a 
constant and v defined by Eq. (|l3|). This rapid decrease 
with increasing temperature is consistent with exact di- 
agonalization studies at hight temperaturesEI. Kliimper 
et al. have found on the other handc3, with an alternative 
Bethe-ansatz approach, a functionally different behavior 
for D(T), see Fig. [|. For a numerical probe of D(T) we 
focus on v = 6 and consider several small temperatures. 
In Fig. (Q) we show a comparism-of our data with the 
two available analytical resultsElro. Our results agree 
with the temperature-dependence predicted by Kliimper 
et al, suggesting in ptarticular a finite Drude weight also 
at the isotropic pointE3. 

Diffusive transport - We consider now a perturba- 
tion of H^ xxz ^ that breaks the integrability but keeps the 
^-component of the magnetization conserved: 



H' = J' Z J2-S? S ; 



i+3 



(14) 



D s is the spin-diffusion constant. Eq. ( |l5| ) is consistent 
with D s = c\ s , where X s = ct is the mean free length. 

We have evaluated <r(T), in addition to the data pre- 
sented in Fig. ||, for v = 6, J z = 0.3 J xx and find a(T = 
0.0047^) = 13.6±0.9, a(T = 0.008J XX ) = 12.1 ±1.0 and 
a(T = 0.012 J xx ) = 10.1 ±0.8. 

We now take J xx — 2000 K, which is appropriate 
for Sr2Cu03 and evaluate the transport coefficients for 
J' z = 0.3 J xx , v = 6, T = 0.004 J xx . We find A w 88 lat- 
tice constants and D s w 6xl0 15 sec -1 . We do not expect 
this model to be directly relevant for Sr 2 Cu0 3 . But if we 
ask ourselves now the question, whether the experimental 
results for A and D s being large (but finite) could be ex- 
plained within a pure spin-model with a small deviation 
from integrability, we might expect H = H^ xxz ^ + H' 
to show the characteristic behavior of a non-integrable 
quantum-spin chain. If we now change the deviation from 
integrability (controlled by J' z ) such that 71 (0) decreases 
by a factor of about ten, then both A and D s would in- 
crease by the same factor and wouldJje consistent with 
the experimentally measured valueaj'H. 

Discussion - We have shown that certain transport 
properties of quantum-spin chains can be evaluated di- 
rectly from two-point correlations functions using a rela- 
tion in between S(q,U) n ) an d the the dyamical conduc- 
tivity which we have derived. For the integrable chains 
we support the original suggestion by Zotos et alH~u of 
a finite Drude weight at finite temperatures and settle 
a recent dispute regarding the functional form of D{T). 
In addition we presents results suggesting the absence 
of ballistic transport (i.e. zero Drude-weight) for a non- 
integrable model, for which we are able to estimate mag- 
nitude of the DC-conductivity. We have shown that the. 
experiments on quasi one-dimensional spin-compoundaH'cl 
are, in principle, consistent with the notion that they 
probe directly non-integrability effects, though we can- 
not rule out at this point that disorder or the coupling 
to phononsQ would lead to the observed finite values for 
the transport coefficients. 

Acknowledgments - This work was supported by the 
DFG. We would like to thank A. Kliimper and X. Zotos 
for giving us their numerical results and their valuable 
comments. We also want to thank A.Sologubenko and R. 
Valenti a careful reading and several suggestions. 



The expression (Q) for the spin current remains valid. 
For this model, H = if' 112 ' + H', we predict a transi- 
tion to a gapped phase around J' z = 0.3 J xx (for v = 6), 
see Fig. [|. We find the relaxation time r = 1/271(0) = 
lim g ^o 1/271(5) to be finite within the numerical accu- 
racy (due to finite-^ and u> n resolution), leading to a finite 
DC-conductivity in the gapless phase. 

For l/u> ^> r the optical conductivity takes (for small 
07/71(0)) the diffusion form 



a(q,uj) 



a uj 



lD s q 2 



27i(0) 



C 2 T. 



(15) 
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FIG. 1. The parameters D^q) and Ai(g) from Eq. (|) 
as a function of momenta q, for the XX-model (J z = 0) and 
various system sizes at T — 0.004-Jjn. Statistical error bars 
are of the order of the system sizecj. 
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FIG. 2. Z>i(g), Ai(g) and 71(g) from Eq. (|) as a function 
of momenta q for the XXZ-model, L — 512 and for various J z 
at T = 0.004 J xx . 71 (q) is too small for J z < J xx to show up 
on this scale. The lines are the Bethe Ansatz result ( |l2| ) for 
the velocity c(J z ) (no fit, for J z < J xx ). For the discussion of 
the fit for J z = 1.5J XX see the text. 
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FIG. 3. For L = 512 and T = 0.004 J xx the QMC-results 
for the gap Ai(0) and the relaxation rate 71 (0) as a func- 
tion of J' z for H {xxz) + H' with J z = J M cos(7r/6). In- 
set: The QMC-results for the Drude weight for J' z — 
and T = 0.004 J xx and L = 512 in comparison with the 
Bethe- Ansatz result at T = 0, Eq. Note that D(T) 

is smooth for low T (see Fig. ^). 
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FIG. 4. QMC results 

L = 512, 1024 and J z = J s 



for the 
M cos(7r/6) 



Drude weight for 
as a function of tem- 



perature (in units of J X x)riQ comparison with two (solid lines: 
Ref.E3, dashed lines: Refo) Bethe-Ansatz results. 



5 



